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5 : Abstract 

X 

I ' Single and multi-photon events with missing energy are analysed using data 

^' collected with the L3 detector at LEP at a centre-of-mass energy of 189 GeV, for 

^ ! a total of 176 pb~^ of integrated luminosity. The cross section of the process e"'"e~ 

-^ uu'-fl'-f) is measured and the number of light neutrino flavours is determined to 
be Ni, = 3.011 ± 0.077 including lower energy data. Upper limits on cross sections 

;J! ' of supersymmetric processes are set and interpretations in super symmetric models 



provide improved limits on the masses of the lightest neutralino and the gravitino. 
Graviton-photon production in low scale gravity models with extra dimensions is 
searched for and limits on the energy scale of the model are set exceeding 1 TeV for 
two extra dimensions. 
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1 Introduction 

In the Standard Model T[| single or multi-photon events with missing energy are produced via 
the reaction e"'"e~ -^ vi"^{'-)) which proceeds through s-channel Z exchange and t-channel W 
exchange. Searches for single and multi-photon final states as well as measurements of the 
e"'"e~ -^ uu'-f{'j) cross section have already been performed by L3 [2H1] and by other LEP 
experiments jS] at lower centre-of-mass energies. For the first time, the determination of the 
number of light neutrino species from single photon events at energies above the Z resonance 
is reported here. 

In supersymmetric models [H] different supersymmetry (SUSY) breaking mechanisms lead 
to different phenomenologies. The SUSY breaking scale, vF, or equivalently the gravitino 
mass (rriQ = F/[y^3/{87r) mp] where mp is the Planck mass), is considered as a free parameter. 
Three different scenarios are distinguished: heavy, light and superlight gravitinos. 

In gravity- mediated SUSY breaking models (SUGRA) the gravitino is heavy (100 GeV < 
itIq < 1 TeV) and thus does not play a role in production or decay processes. The lightest 
neutralino is the lightest supersymmetric particle (LSP), which is stable under the assumption 
of R-parity [Zj conservation and escapes detection due to its weakly interacting nature. In 
this scenario, single or multi-photon signatures arise from pair-production of neutralinos (X1X2 
and X2X2) 0. Subsequent one-loop decays of X2 into X17 have a branching fraction close to 
100% if one of the two neutralinos is pure photino and the other pure higgsino |9|. In general, 
neutralinos are mixtures of photinos, zinos and higgsinos. 

In models with Gauge-Mediated SUSY Breaking (GMSB) [101, a light gravitino (10"^ eV < 
777.(5 ~ -'-0^ ^) ^^ ^^^ LSP. In this case the gravitino plays a fundamental role in the decay of 
SUSY particles. In particular, the Xi is no longer stable and decays through Xi — > G7 if it 
is the next-to lightest supersymmetric particle (NLSP) [11^. Pair-production of the lightest 
neutralino leads to a two-photon plus missing energy signature in the detector. 

When the scale of local supersymmetry breaking is decoupled from the breaking of global su- 
persymmetry as in no-scale supergravity models ^^, the gravitino can be superlight (10~^ eV < 
rriQ < 10""^ eV). Then, it is produced not only in SUSY particle decays but also directly in 
pairs ^H! or associated with a neutralino ^1]. Pair-production of gravitinos accompanied by 
initial state radiation leads to a single photon signature. This signature also arises in XiG 
production when the neutralino decays radiatively to gravitino and photon. 

Recently, it has been proposed that the fundamental gravitational scale in quantum gravity 
models with extra dimensions is as low as the electroweak scale fSj thus naturally solving the 
hierarchy problem. Within the framework of these models real gravitons are produced in e"'"e~ 
collisions through the process e"'"e~— »-7G, where the graviton escapes undetected leading to a 
single photon plus missing energy signature. 

2 Data Sample and Simulation 

In this analysis we use the data collected by the L3 detector TB^ during the high energy run of 
LEP in 1998 corresponding to an integrated luminosity of 176.4 pb~^ at an average centre-of- 
mass energy of ^/s = 188.6 GeV, hereafter denoted 189 GeV. 

Monte Carlo events for the following Standard Model processes are simulated: e^e~ — >■ 
i/z/7(7) with KORALZ [T7], e"'"e~ — *> 77(7) with GGG [IH], Bhabha scattering for large scattering 
angles with BHWIDE ^\, and for small scattering angles with TEEGG [201; and four-fermion final 



states specifically the processes e"'"e —>■ e'^e e'^e with DIAG36 J2T], and e"'"e —>■ e'^e uv with 
EXCALIBUR [22]. 

SUSY processes are simulated with the Monte Carlo program SUSYGEN |2S1 for SUSY particle 
masses (ttisusy) between zero and the kinematic limit and, in Xi LSP scenarios, for Am = 
"^susY — '"^LSP between 1 GeV and msusY- To ensure the radiative decay of the neutralino, 
the scalar electron (cr) mass is set to 100 GeV, except for XiG production where it is set to 



200 GeV. The detector response is simulated using the GRANT program |23], which takes into 
account the effects of energy loss, multiple scattering and showering in the detector. 

3 Event Selection 

Electrons and photons are measured accurately by the BGO electromagnetic calorimeter. They 
are required to have an energy greater than 0.9 GeV. The shape of their energy deposition 
must be consistent with an electromagnetic shower. Electrons are defined as electromagnetic 
clusters matched with a charged track reconstructed in the central tracking chamber. Identified 
conversion electrons from photons interacting with the beam pipe or with the silicon microvertex 
detector, 4% of the total, are also accepted as photons. Bhabha events and e"'"e~ —>■ 77(7) events 
that are fully contained in the calorimeter are used to check the particle identification as well 
as the energy resolution, which is 1% for high energy electrons and photons in both the barrel 
and the endcaps. The barrel region is defined as the polar angle range 43° < 9 < 137° with 
respect to the beam axis and the endcap region as the polar angle range 14° < 6 < 36° or 
144° <e < 166°. 

3.1 High Energy Photons 

The selection of high energy single and multi-photon events requires at least one photon with 
energy greater than 5 GeV in the barrel or endcaps region. There must be no charged tracks 
apart from those consistent with photon conversion. The following cuts are imposed to suppress 
events which do not consist of photons only in the final state. The energy not assigned to iden- 
tified photons has to be smaller than 10 GeV and the energy measured in the electromagnetic 
calorimeter between BGO barrel and endcaps must be smaller than 7 GeV. There must be no 
track in the muon chambers and at most one BGO cluster not identified as a photon. 

To reduce the background from radiative Bhabha events with particles escaping along the 
beam pipe, as well as from the process e"'"e~ —>■ 77(7), events with less than 20 GeV transverse 
momentum are rejected if energy is observed in the small polar angle detectors covering an 
angular range of 1.5° — 11°. The total transverse momentum of photons is required to be greater 
than 5 GeV if no second photon with energy greater than 5 GeV is found. If two calorimetric 
clusters are present and if only one is identified as a photon, their acollinearity must be greater 
than 5.2° and their acoplanarity must be greater than 2.4°. Furthermore, energy clusters in 
the hadron calorimeter (HCAL) must have less than 3 GeV energy if a photon is detected with 
an acoplanarity less than 15° to the HCAL cluster. 

When a second photon with energy above 5 GeV is present, the total transverse momentum 
must be greater than 3 GeV and the recoil mass must be larger than 20 GeV. If the total 
transverse momentum is smaller than 30 GeV, the acollinearity is required to be larger than 
8.1° and the acoplanarity to be larger than 5.2°. If the transverse momentum is smaller than 20 
GeV, the missing momentum direction is required to be at least 7° away from the beam pipe. 
If the acoplanarity is smaller than 2.4°, the recoil mass must be greater than 50 GeV. 



To suppress cosmic ray background, we require for photon energies smaller than 15 GeV, 
that the most energetic photon is not aligned with hits in the muon detector. For photon 
energies larger than 15 GeV, there must be at least one scintillator time measurement within 
±5 ns of the beam crossing time. Furthermore, an event is rejected if more than 20 hits are 
found in the central tracking chamber in a 1 cm road between any pair of energy depositions 
in the BOO. 

The number of events with one or more photons is listed in table Q together with the 
predicted rates for vi"-){'-)) and other processes originating from e"'"e~ collisions. The cosmic ray 
background in the event sample is estimated from studies of out-of-time events and also listed 
in tabled Figure ^ shows the energy spectrum of the most energetic photon normalised to the 
beam energy for single and multi-photon events. 

For the sub-sample of events with two or more photons a minimum energy for the second 
photon of 1 GeV is required. In the data 21 events are observed compared to a Monte Carlo 
prediction of 36.2 events, see table ^ For recoil masses larger than 110 GeV we observe 2 
events compared to an expectation of 12.7 events. Figure E^a) shows the two-photon recoil 
mass distribution. The lack of data compared to the Monte Carlo prediction has been subject 
to extensive investigations concerning the performance of sub-detectors and triggers used in 
this analysis. The noise level of sub-detectors is studied using randomly triggered beam-gate 
events. The performance of the electromagnetic calorimeter is cross-checked with Bhabha 
events and events from e"'"e~ — >■ 77(7). Triggers important for single and multi-photon events 
are investigated using single electroio and Bhabha events. The theoretical predictions for the 
cross section of e'''e~ — > uu'-f'-fl'-f) obtained from KORALZ and NUNUGPV are found to agree within 
5%. No systematic effect is found to explain the low two-photon rate. An independent analysis 
leads to the same conclusion. Therefore the deficit is treated as a statistical fluctuation. 

3.2 Low Energy Photons 

This selection extends the energy range for photons down to 1.3 GeV. It covers only the barrel 
region where a single photon trigger is implemented with a threshold around 900 MeV |25j . 
There must be no other BGO clusters in barrel or endcaps with more than 200 MeV. The energy 
in the HCAL must be less than 6 GeV. To reduce the rate of small angle Bhabha scattering no 
energy deposit is allowed in the forward detectors. Events with a track in the central tracking 
chamber or in the muon chambers are rejected to reduce the rate of single electron and cosmic 
muon events. To further reduce cosmic ray events not pointing to the interaction region, cuts 
on the transverse shape of the photon shower are made. For the simulation of the process e"'"e~ 
— >■ e"'"e~7(7) the TEEGG program is used which includes fourth order contributions. In order to 
use this program, a cut on the transverse momentum of the photon greater than 1.3 GeV is 
applied [^. From a study of single electron events its precision is estimated to be at the 20% 
level. The number of selected events, predictions from e"'"e~ collision processes and an estimate 
of cosmic ray contamination are listed in tabled Figure |21^b) shows the observed photon energy 
spectrum compared to Monte Carlo prediction. 



^'Radiative Bhabha scattering events where one electron and a photon have a very low polar angle, and only 
a low energetic electron is scattered with a large polar angle. 



4 Neutrino Production 

To measure the cross section of the e"'"e~ — >■ i>u'j{'y) process we restrict the analysis to photon 
energies above 5 GeV (see table H} to ensure a good signal to background ratio. 

The overall efficiency for e"'"e~ — > i>u'y{'y) events satisfying the kinematic requirements E^ > 
5 GeV and |cos^^| < 0.97 is 60.77%. This efficiency includes a correction of (2.58 ± 0.18)% 
due to cosmic ray veto requirements, which is estimated by studying single electron events. A 
correction of (0.67 ± 0.07)% due to detector noise sources not properly simulated such as that 
induced by beam halo in the forward detectors is also included and is quantified using randomly 
triggered beam-gate events. The systematic error on the efficiency is composed of the errors 
on the two corrections and several other sources including an error of 0.34% due to photon 
identification, of 0.60% due to an uncertainty on the amount of converted photons, and 0.21% 
due to limited Monte Carlo statistics. The total systematic error on the efficiency amounts to 
0.75%. The error on the luminosity is 0.4 pb~^ and on the total background contamination the 
error is 2.1 events. 

The measured cross section at y/s = 188.6 GeV is 

c^up-yi-y) = 5.25 ± 0.22 (stat) ± 0.07 (syst) pb 

to be compared to the prediction of the Standard Model of 5.28 ± 0.05 pb obtained with KORALZ, 
where the 1% error accounts for the theoretical uncertainty assigned to this process [221 ■ This 
measurement is extrapolated to a total cross section for e"'"e~ — >■ vi'{'^) production of 58.3 ± 
2.7 pb. The prediction of the Standard Model obtained with KORALZ is 58.6 pb. Figure 01 shows 
both the vi"-){'-)) cross section measurement and the total neutrino-pair extrapolation versus 
centre-of-mass energy together with the prediction of the Standard Model and measurements 
at lower centre-of-mass energies. 

To determine the number of light neutrino species a maximum likelihood ffi to the photon 
energy spectra is performed at each centre-of-mass energy above the Z resonance. For each 
energy interval the theoretical prediction is obtained by linearly interpolating KORALZ predic- 
tions for Ny = 1,2,3,4,5. Due to the different contributions to the energy spectrum from 
UeUe t-channel production via W exchange and ui? s-channel production via Z exchange, this 
method is more powerful than using the total cross section measurement. In addition to the 
systematic error from the cross section measurement, the theoretical uncertainty on the photon 
energy spectrum - estimated by comparing KORALZ with NUNUGPV - is taken into account. The 
result is N^ = 3.05 ± 0.11 ± 0.04. 

A compilation of the measurements at the different centre-of-mass energies is shown in table 
121 The precision of this result is comparable with our previous measurement ^ from single 
photon events around the Z resonance. The combined measurement is 

A/'^ = 3.011 ±0.077. 

This result is more precise than the present world average on the number of light neutrino 
families determined with the single photon method |27j . 

5 Limits on Super symmetry 

The limits derived in the following are obtained from 189 GeV data. They correspond to a 
confidence level of 95%. Cross section limits are calculated using a likelihood approach 



where the spectra of a discriminant variable for data, background and signal simulations are 
compared. 

5.1 Single Photon Signature 

For interpretations within SUSY frameworks the low energy photons are included. Here, the 
discriminant variable used to derive cross section limits is the photon energy. 

In the heavy gravitino scenario, the single photon signature arises from the reaction e"'"e~ 
-^ X2X1, which proceeds through s-channel Z exchange and t-channel scalar electron exchange 
(gl.r)- Cross section upper limits shown in figure IH are set under the assumption of 100% 
branching fraction for X2 — > X17. Typical detection efficiencies for this process are around 75%. 

Also the reaction e"'"e~ -^ GXi proceeds through s-channel Z exchange and t-channel eL,R 
exchange. Efficiencies for this process with Xi -^ G7 range between 64% for m^o =0.5 GeV 
and 79% for m^o at the kinematic limit. The cross section upper limit as a function of m^o 
is shown in figure El^a) together with the expected limit obtained in Monte Carlo trials with 
background only. The no-scale SUGRA model of fl], referred to as LNZ, has only two free 
parameters - gravitino and neutralino masses. The neutralino is the NLSP, which is almost 
pure bino. Here, the dominant decay channel is Xi -^ G7. The small contribution of the decay 
into Z for m^o > 100 GeV is taken into account. Figure Efb) shows exclusion contours in the 
rriQ-m^o plane. 

If mNLSp > y/s, the process e"'"e~ — *> GG [13j is the only reaction to produce SUSY particles. 
Accompanied by initial state radiation it leads to single or multi-photon signatures. Following 
our analysis in |3j a lower limit on the gravitino mass is derived 

m^ >8.9- 10"^ eV, 

corresponding to a lower limit on the SUSY breaking scale of a/F > 192.3 GeV. The average 
lower limit for the gravitino mass obtained in Monte Carlo trials with background only is 
9.7 ■ 10-6 eV. 

5.2 Multi-Photon Signature 

Using a binned likelihood technique, the discriminant variable is constructed for the multi- 
photon events combining the energies of the two most energetic photons, their polar angles, 
recoil mass, and the polar angle of the missing momentum vector. For each event class, back- 
ground and signal Monte Carlo processes, denoted by j, and each input quantity i, a probability 
density function /* is computed and the discriminant variable of an event is then given by 



F{x) = ^";^^ , \ with p){x,) 



where Xi are the measured values of the six input variables of an event. The distribution of 
the discriminant is shown in figure IHl for XiXi — > GG77 with m^o = 90 GeV. The discrepancy 
between measurement and Standard Model prediction is located in the background and not in 
the signal region. This holds also for the other mass points and for the heavy gravitino scenario. 
In the heavy gravitino scenario, a two-photon signature is produced by the process e'^'e" — *■ 
X2X2 and subsequent decay to X17. Typical efficiencies for this process are around 64%. Cross 
section upper limits are obtained as shown in figure [TU^a). The interpretation of the ee77 event 



with large transverse missing energy observed by CDF 122] suggests a high branching ratio for 
the radiative decay of the X2 in the heavy gravitino scenario, which can be achieved if X2 is a 
pure photino and Xi is a pure higgsino. With this assumption, the lower mass limit of X2 as a 
function of the scalar electron mass is calculated for mass differences between X2 and Xi greater 
than 10 GeV as shown in figure El^b). For each X2 mass, the exclusion is obtained using the 
most conservative cross section upper limit for any Am > 10 GeV. The regions kinematically 
allowed [^ for the CDF event are overlayed in figure [7|^b). The two exclusions obtained for 
equal masses of eL,R and for rrie^ 3> rrie^ are shown in the interesting mass range for rrie^. 

The selection described in this paper is devised for photons originating from the interaction 
point. For a neutralino mean decay length d^^o larger than 1 cm the experimental sensitivity 
drops. This problem arises only for peculiar situations in the light gravitino scenario. The 
following limits in the gravitino LSP scenario are derived under the assumption of d^o < 1 cm. 

The cross section limits for the process e"'"e~ -^ XiXi -^ GG77 are displayed in figure IHl^a) 
versus neutralino mass. The efficiency varies between 17% for m^o = 0.5 GeV and 62% for 
m^o = 94 GeV. Theoretical predictions for two extreme cases of neutralino contenio |31j . 
which determines its coupling to the photon, are shown in the same figure. For these cases of 
neutralino composition and for a pure photino, we derive lower limits on the mass of the lightest 
neutralino as listed in table El Figure Efb) shows the exclusion in the Xi — eL,R mass plane 
derived with our data for a neutralino being pure bino. The ee77 event observed by CDF also 
has an interpretation in supersymmetric models with gravitino LSP [21]. Our analysis almost 
rules out this interpretation as shown in figure IHfb). 

In minimal models with gauge-mediated SUSY breaking only five parameters determine 
the sparticle sector of the theory jTU]. The parameters are A, the scale of SUSY breaking 
in the messenger sector, M^, the messenger mass scale, Nm, the number of messenger fields, 
tan/3, the ratio of Higgs vacuum expectation values. In this model the absolute value of /i, 
the higgsino mass term, is fixed, however its sign is a free parameter. They have been scanned 
to obtain neutralino masses, pair-production cross sections and branching ratios for the decay 
to gravitino and photon. The scan ranges on the individual parameters are [S2] 10 TeV < 
A < 100 TeV, A/0.9 < M„ < A/0.01, A^;^ = 1...4, 1 < tan/3 < 60, sign/i = ±1. The 
program ISASUSY [^ has been used to calculate sparticle masses and couplings from GMSB 
model parameters, and SUSYGEN to derive from these numbers the cross section for neutralino 
pair-production including initial state radiation. Assuming a neutralino NLSP scenario, the 
minimal cross section of XiXi production obtained within GMSB is shown in figure |H^a), which 
leads to a lower limit of 

m^o > 88.2 GeV . 

6 Limits on Graviton Production 

Massive spin 2 gravitons propagating in 4+6 dimensions interact with Standard Model particles 
with sizable strength in low scale gravity models with extra dimensions J3] . Gravitons produced 
via e"'"e~ -^ 7G lead to a single photon and missing energy signature, since the graviton is not 
observed in the detector. The reaction proceeds through s-channel photon exchange, t-channel 
electron exchange and four-particle contact interaction [HI]. 

^^For the higgsino case a 2% photino component is required to ensure the decay into 7G. 



To convert the theoretical cross section of this process |21| into an estimate on the number 
of events expected from graviton production, the differential cross section in energy and angle 
has been multiplied by efficiency and luminosity. The efficiency is derived from vi"^{'^) Monte 
Carlo simulation in a grid in the E^ — cos 9^ plane for E^ > 4 GeV. The efficiency for e"'"e~ — *■ 
7G within E^ > 4 GeV and cos 6'^ < 0.97 is listed in tablelUfor 2 < 5 < 10. The energy spectra 
shown in figure Wljo) and El^b) are used to derive upper limits on the cross section. They are 
listed in table IH together with the corresponding values for the energy scale Mc- These bounds 
improve on our previously published limits |SS]- 
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E^>h GeV 


E^ > 1.3 GeV 
E^<5 GeV 


E^, > 5 GeV 
E^2 > 1 GeV 


Total 


Barrel 


Endcaps 


Data 


572 


297 


275 


395 


21 


z/z77(7) 


567.3 


288.9 


278.4 


48.7 


35.5 


e+e" background 


6.5 


2.2 


4.3 


358.5 


0.7 


Cosmic background 


3.1 


1.1 


2.0 


3.6 





Total expectation 


576.9 


292.2 


284.7 


410.8 


36.2 



Table 1: Number of events selected in data, Monte Carlo predictions for processes 
from e"'"e~ collisions and contamination of cosmic ray background in the indicated 
kinematic regions. 



v^ (GeV) 


N, 


130.1 


2.63 ± 0.40 ± 0.10 


136.1 


2.98 ± 0.49 ± 0.14 


161.3 


3.68 ± 0.53 ± 0.09 


172.1 


4.24 ± 0.65 ± 0.09 


182.7 


3.13 ± 0.26 ± 0.05 


188.6 


2.94 ± 0.15 ± 0.04 


130 - 189 


3.05 ± 0.11 ± 0.04 


88-94 


2.98 ± 0.07 ± 0.07 


Average 


3.011 ± 0.077 



Table 2: Number of neutrino families measured from single photon events. 



Xi content 


"^eL.R (GeV) 


ml^ (GeV) 


Bino 


150 


87.9 


Bino 


100 


90.8 


Photino 


150 


88.3 


Photino 


100 


91.1 


Higgsino 


— 


89.0 



Table 3: Neutralino mass limits for several neutralino compositions and selectron 
masses. 



5 


2 


3 


4 


5 


6 


7 


8 


9 


10 


e(%) 


42.8 


40.7 


38.9 


37.6 


36.5 


35.5 


34.7 


34.0 


33.4 


<^ (pb) 


0.638 


0.646 


0.651 


0.658 


0.664 


0.670 


0.674 


0.678 


0.680 


Md (GeV) 


1018 


812 


674 


577 


506 


453 


411 


377 


349 



Table 4: Selection efficiency e for e"'"e — ^ 7G, upper cross section limit and lower 
limit on the energy scale Md as a function of the number of extra dimensions S. 
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Figure 1: (a) Energy of the highest energetic photon normahsed to the beam en- 
ergy for single and multi-photon events at 189 GeV in the barrel region, (b) Same 
distribution with the endcaps included. 
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Figure 2: (a) Recoil mass for the multi-photon sample, (b) Low energy part of the 
energy spectrum of single photon events. 
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Figure 3: Production cross section of e"'"e~ -^ ^^(l) and e^e~ — > z/z/7(7) as a 
function of the centre-of-mass energy. Points with error bars represent the z/z/7(7) 
measurements and squares with error bars are the extrapolation to I'i^i'y). The full 
line is the theoretical prediction for N^, = 3 and dashed lines are predictions for 
A^,, = 2, 4 as indicated. 
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Figure 4: Upper limits on the production cross section in picobarn for the process 
e+e~ — > X2X1 -^ X1X17 assuming 100% branching ratio for X2 -^ X17. 
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Figure 5: (a) Upper limits on the production cross section for the process e+e^ -^ 
GXi — »• GG7 and average hmit obtained using Monte Carlo trials with background 
only, (b) Region excluded in the LNZ model in the plane uIq versus m^o. 
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Figure 6: Discriminant variable for XiXi -^ GG77 with m^o = 90 GeV. The signal 
corresponds to the upper limit of 3.15 events derived for this mass point. 
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Figure 7: (a) Upper limits on the production cross section in picobarn for the process 
e+e~ — > X2X2 — > X1X177. (b) Excluded region in the neutralino selectron mass plane. 
The shaded region corresponds to mg^ ^ rrie^ and the hatched region is additionally 
excluded when mgL = '"^ 
as a function of m^o are indicated, where Xi 



Regions kinematically allowed for the CDF event [SO] 
Hfe = H?sin/5 + HOcos/3. 
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Figure 8: (a) Upper limit on the cross section for e"''e 
predictions for two extreme cases of X\ composition and for the most conservative 
GMSB prediction are also shown, (b) Excluded region for a pure bino neutralino 
model compared to the region consistent with the supersymmetric interpretation of 
the CDF event in the scalar electron scenario |31j . 
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